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RADIOFREQUENCY POWER TRANSFER TO ION-CYCLOTRON WAVES 


IN A COLLISION-FREE MAGNETOPLASMA 
by Donald R. Sigman 
Lewis Research Center 


SUMMARY 

With the Fourier integral theory developed by T. H. Stix, calculations 
were made on power transfer from a radiofrequency coil to an ion-cyclotron wave 
in a cold collisionless cylindrical plasma. The effects of varying parameters 
such as static magnetic field, electron density, radiofrequency, plasma and coil 
radius, and wavelength are discussed. It is shown that there is an optimum 
radiofrequency coil wavelength for maximum power transfer, and curves are pre- 
sented to facilitate the selection of this wavelength. The importance of 
placing the radiofrequency coil as close to the plasma as possible is also 
shown. Finally, a comparison is made between calculations from the above men- 
tioned Fourier integral theory and a former Fourier-Bessel series theory (also* 
developed by Stix). The two theories are not found to be in agreement. 


INTRODUCTION 

The generation of ion-cyclotron waves in a plasma column is of interest to 
experimentalists concerned with the heating of plasma through collisionless 
processes. This report contains calculations from a theoretical plasma model 
of the effect of various plasma parameters on the power transferred from a 
radiofrequency coil to the ion-cyclotron wave in the plasma. In this model, the 
radiofrequency coil is of finite length, and it is wrapped around an infinitely 
long plasma column surrounded by a vacuum and immersed in a steady uniform 
magnetic field parallel to the axis of the column. The assumption that the 
plasma is of constant density, cold, and collision free is made in the model. 

In the region underneath the coil, power is transferred from the radio- 
frequency coil to the ion-cyclotron wave. This wave then propagates unattenu- 
ated out both ends of the coil and down the plasma column. Since the plasma 
model is cold, all energy from the coil shows up as electromagnetic energy of 
the wave. 

By using the previous physical model and assuming, in addition, that the 
electron has zero mass, T. H. Stix (ref. l) derived expressions for the power 
transferred to ion-cyclotron waves. 



Stix, by using a Fourier integral, represented the current sheet distribu- 
tion in the finite length coil by an infinite series of current sheets, each 
infinitely long. The wavelengths of these current sheets ranged from zero to 
infinity, and the relative strength of each current sheet was determined by the 
Fourier integral. Since the electric fields associated with each of these in- 
finitely long current sheets can be easily calculated, it was then possible to 
get an expression for the field of the finite length coil by superposing the 
fields of all the infinitely long current sheets. This technique resulted in 
an integral for the electric field of the finite length coil. This integral 
was solved by using contour integration to arrive at the final expressions for 
both the real and imaginary parts of the field (i.e., those parts of the field 
in phase with and 90° out of phase with the current in the current sheet). The 
power transferred to the wave was then computed by averaging over the length of 
the coil the product of the current in the finite current sheet and the real 
part of the electric field at the current sheet. 

Stix formerly used another technique that employed a Fourier series to cal- 
culate power transfer to the ion-cyclotron wave (ref. 2). A comparison of the 
results of this former technique with those of the .Fourier integral technique 
is given in the appendix. 

In this report, the results of digital computer calculations from the ex- 
pressions derived in reference 1 are presented. The effects of the following 
plasma parameters on power transfer were studied: plasma wavelength, plasma 

mode, radiofrequency coil radius and length, and magnetic field. An attempt 
has been made to present the results of the calculations in such a way that 
they are an aid to an experimentalist both in designing an efficient means for 
power transfer from a radiofrequency coil to the plasma wave and in interpret- 
ing the data. 


SYMBOLS 

The cgs-gaussian units are used throughout, 
a length of one-half of radiofrequency exciting coil 

B magnetic field strength 

Bq magnetic field at ion-cyclotron particle resonance 

c velocity of light 

Eg see eq. (4b) 

eY 6 component of vacuum electric field, see eq. (4a) 

U 

e electronic charge 

F radiofrequency coil radius correction factor 

G defined by eq, (12b) 
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Iq modified Bessel function of first kind 

Jq Bessel function of first kind 

magnitude of radiofrequency coil current 
j radiofrequency coil current per unit length 

Kq modified Bessel function of second kind 

k axial wave propagation constant (wave no,) in plasma 

k n see eqs. (8) and (9) 

koi axial radiofrequency coil current propagation constant (wave no.) 
mj mass of hydrogen ion 

n e electron density 

P^ total power transferred to ion-cyclotron waves from radiofrequency coil 

P* relative power transfer (2jta 2 j^P* = P^ ) 

p plasma radius 

R relative power transfer factor, see eq. (6d) 

r radial cylindrical coordinate 

S^(rj) shape factor for imaginary power transfer, see eq. (6a) 

S r (ri) shape factor for real power transfer, see eq. (6b) 
s radiofrequency coil radius 

t time 

z axial cylindrical coordinate 

T) 2a(k oz - k n ) 

0 azimuthal cylindrical coordinate 

v radial wave number 

Hq ion plasma frequency 

co radiofrequency/2 it 
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THEORY 

The axis of the plasma cylinder is designated the z-axis and cylindrical 
coordinates are used. The steady magnetic field is directed along this z-axis. 
The exciting radiofrequency coil for the system has a length 2a, surrounds the 
plasma cylinder, and is centered at z = 0 (see fig. l). The aximuthal cur- 
rent per unit length in the coil is to be represented by the following Fourier 
series : 




1=1 


i^Q^z-cot ) 


-a < z < a 


(la) 


j* ■ 0 | z | > a (lb) 

In order to represent a finite length coil in a single equation, Stix fur- 
ther decomposed each Fourier current component in equation (l) into its Fourier 
integral spectrum so that 


m 

3* = J2 
1=1 
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This led to the following expression for the electric field Eg at the radius 
of the coil: 


where 


m 

(®@) = S — j j 

° r=s / 4 rt b 

1=1 

E 0 r=s = fsIpCksjK^ks) + Eg 

(A- 

' ' c 2 kJ ^('vp )K^( kp ) - vj|(vp)K 1 (kp) 


/*00 

/ (C, 

-00 


sin(k - 


k 0Z 


(k - k Q £ ) 


e i(kz-cst) di 


(3) 


(4a) 

(4b) 


The 




is the 9 component of the vacuum electric field at r = s when 


the coil is assumed to be infinitely long and to have a current distribution 
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per unit length j ^ e 1 (kz-a)t) . gy working from equation (3), Stix was able to 
derive the following expression for the power transferred to the n kil natural 
mode of the plasma wave when all the j j 1 s have been assumed zero except one: 


-i(2*sa2)jj[(k - >%)^ r , s ] k=j£n [s r (,) - iSjU)] 


+(2na 2 )jfE[s r (i,) - lS ± ( n )] 

(5) 

S r (r) ) = 2 .(.. 1 ~ cos t ll 
i 2 

(6a) 

s ± (-n) = ha - jiB. al 
*1 

(6b) 
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(6c) 

R = -i(k - k,.) (eI ) — 

^ V e r =s/k=k n 

( 6d) 


The real part of P^ corresponds to resistive plasma loading, which is the 
power transferred to the plasma wave as measured in the laboratory. Thus, 

Re P z = 2jra 2 jfRS r ( 11 ) ( 7 ) 

The natural modes are determined by finding the conditions under which Eg 
in the plasma under the infinitely long coil remains nonzero when j* = 0. This 
condition is met when the denominator of equation (4b) is zero: 

kJ x ( vp)Kj_(kp) - vJ£( vp)K 1 (kp) =0 (8) 

is given by 


( 9 ) 


The simultaneous solution of equations (8) and (9) gives the values v n 
and k n for the n^* 1 natural mode of the plasma. These values were then used 
in equation (6d) to compute R on the 7094 digital computer. At k = k^, both 
numerator and denominator of equation (6d) are zero. This in determ ina ncy is 
eliminated in the calculations through use of 1 'Hospital's rule with respect 
to k. 


The dispersion relation for ion-cyclotron waves in the plasma 
equations (5) to (9) in reference 1: 
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DISCUSSION OF RESULTS 


For the purpose of the discussion that follows, it should be noted that 
the equation for power transfer (eq. (7)) contains three terms of interest: 
2«a 2 j^, R, and S r (r]) . The function S r (r]), which is called the "shape factor" 
by Stix, is the only term which contains k()7> an d it is a maximum when 
kQ i = k n (see fig. 2). The term R is a function of k n , s, p, n e , and G, 
but not kQj . Thus, the maximum RS r (r}) product is achieved by selecting k qi 
equal to the value of k n for which R is a maximum, since under this condi- 
tion both R and S r (T}).are maximum. 

In the following sections, the effects of static magnetic field, wavelength 
and radii of the plasma and radiofrequency coil, and plasma density on radio- 
frequency power transfer are discussed. Also, plasma modes other than the 
primary mode are compared to the primary mode, and calculations are made with 
nonsinusoidal radiofrequency coil current distributions as well as with sinus- 
oidal distributions. In equation (l), it is stated that the current sheet is a 
traveling wave current sheet in the +z direction. The results presented in the 
following sections are the same if the current sheet is assumed to travel in 
the -z direction. Thus, if one wishes results for a standing wave current 
sheet, he may to a good approximation simply add the contributions from the 
positive and the negative z-directed waves. In doing this, the cross terms, 
which are very small, are neglected. 


Different Plasma Modes 

The simultaneous solution of equations (8) and (9) gives an infinite 
number of values for v n , and for each value of v n there is a contribution to 
power absorption. However, it was found numerically that only the contribution 
for the mode with the lowest value of v n was of significance. In all cases 
studied, the contribution of the lowest mode was more than ten times that of 
the next highest mode; the shapes of the curves presented are not affected by 
inclusion of the higher modes. Thus, all curves shown in this report are from 
lowest order mode calculations (i.e., n = l). 


Magnetic Field 

The dependence of R on the magnetic field for a 0 typical laboratory 
plasma is shown in figure 3. Relative power absorption (P* = RS r (r))) against 
Q is given in figure 4. In this situation, jq = 1 (all other j i = 0), 
kQi = 0.140 per centimeter, a = 2:rt/koq, n e = 5 X 10 H per cubic centimeter, 
plasma radius =2.5 centimeters, coil radius =7.5 centimeters, and 
cd = 2itX6 . 5X10® per second (Bq = 4200 G). The relative power absorption curve 
in figure 4 is very much like a plot of S r (rj) against fi (also presented in 
fig. 4). However, the relative power absorption curve has a narrower half- 
width, and the magnetic field at maximum power absorption is less than the 
magnetic field at the peak in the S r (r}) against G plot (i.e., where r] = 0). 

Even though the calculations for figures 3 and 4 have been made by using 
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specific values of the parameters cp, s, p, kQi, and n e , the basic shape of 
these curves is preserved for other values of the parameters. Different values 
of o, s, or p affect the magnitude of R for a given fl; n e largely deter- 
mines the value of 0, where R is a maximum ( ft decreases as n e increases). 
Different values of kQi move the peak of S r ( t) ) with respect to Q. 

In the derivation for R (eq. (5d) ), the displacement current term from 
Maxwell's equations was neglected. For large magnetic fields (when the phase 
velocity oo/k is no longer negligibly small compared to the velocity of light), 
displacement current may be included by replacing k in the expression for R 

by k 2 - fo z /' z )i /2 . For Bo = 75 kilogauss (co « 2jtX1.16x10 S see”'**) and 
k = 0.15 per centimeter, (co/k ) 2 = 0.24X10 2 ^ = 0.26 c 2 ; therefore, the curve P* 
against Q, is not changed from that of figure 4 where Bq = 4.2 kilogauss. 


Plasma Wavelength and Radiofrequency Coil Length 

In studying the effects of coil length on plasma power absorption, two 
cases are considered,. For case A, the physical length 2a of the coil is held 
constant, and the value of kQp in the principal Fourier component of the cur- 
rent distribution as given by equation (la) is varied. In case B, the physical 
length of the coil is varied but phased for two wavelengths so that a = 27r/k 01 . 

Figures 5(a) to (d) show R against k n for case A for a range of the 
variables p and n e with s held constant. From these curves, it is seen 
that there is a value of k n which makes R a maximum. This point of maximum 
R, however, occurs at different values of k n for different values of p, n e , 
and s. For larger values of s, the peak was found to occur at smaller values 
of k n such that when n e is held constant the product k n s at the peak is a 
constant. The position of the peak is less sensitive to changes in p. For 
high densities (n e > 10^- 2 cm" 3 ), the peak occurs at larger values of k^, but 
it does not vary significantly with density for n e < 10^-2 p e r cubic centimeter. 

2 2 

For case B, Pn = 2«a j£RS r (r[) and a = 2it/koi* When , k 01 = k n ( i . e . , when 
S r ( r] ) is a maximum), then P^ ~ a 2 R = 4jt 2 /k^ R. Figures 6(a) to (d) show 
R/k 2 as a function of for a range of the variables p and n e with s 

again set equal to 10 centimeters. Again there are maxima in the curves, which 
indicate a point of maximum power transfer. The dividing of R by k 2 moves 
the peak to lower k n values (compare figs. 6(a) and (b)). Qualitatively, the 
effects of s, p, and n e on the maximum power points are the same for case B 
as for case A. 

Since figures 5 and 6 are for a fixed value of s = 10 centimeters, a cor- 
rection factor is needed if the curves are to be useful for other values of s. 
We therefore define a function F(k n s) = [k n sK/p(k n s)]2, which contains the 
entire dependence of power on s . This function is plotted in figure 7 for a 
range of s values. Corrections to the values in figures 5 and 6 for values 
of s f 10 may be made by multiplying by the correction factor F(k n s)/F(l0k n ) . 

All the curves in figures 5 and 6 are plotted with k n as the abscissa. 
This is a useful form for designing an experiment. For interpreting data. 
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however, a more useful variable is SI, which is defined by 



1 + 


£_ 


(io) 


Equation (10) permits the evaluation of k n- in terms of measurable quantities 
irj and Si. 


Plasma Density 

Figure 8 shows the effects of varying the density on power adsorption. 

For curve (a), the value of kgi was set at a fixed value (0.140 cm~l) and 
"peak" relative power transfer (P* at k n = kp) was plotted. For curve (b), 
the value of koq was adjusted for each density so that koq was equal to the 
value of k n which gave a maximum in the R against k n curve (fig. 5(a)). 
For curve (c), the value of koi was adjusted to the value of k n which gave 
a maximum in the R/k 2 against k^ curve (fig. 6(a)). 

The conclusion is that over the density range 3x10® per cubic centimeter 
< n e < 5X1012 per cubic centimeter power absorption is a weak function of den- 
sity , but it drops off rapidly if an attempt is made to go to higher densities. 

Figure 9 shows the values of Si for "peak" relative power absorption 
against density when koi = 0.140 per centimeter. 


Plasma and Radiofrequency Coil Radius 

Figure 10 shows the effect of plasma and coil radius on "peak" relative 
power absorption by the plasma. It is seen that when the radiofrequency coil 
radius s is held constant, power absorption rises faster than (p/s) 4 . This 
emphasizes the importance of operating with the radius ratio p/s as near to 
1 as physically possible. 


Other Fourier Current Components 

Because the current distribution of a real laboratory coil cannot realis- 
tically be represented as a single sine wave, it is of interest to compute the 
contribution to power absorption from Fourier components of the coil current 
other than the principal component. For an illustrative calculation it is 
assumed that 

j* = j (o.l51 e i4k 01 z + 0.422 e i2k 01 z + 0.270 e ik 01 z j | z| < a (ll a ) 

j* = 0 z > a (lib) 

where a = 2jt/2koi, koi = 0.070 per centimeter, ne — 5X10^ per cubic centimeter. 
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s/p = 3 and oo = 2 jtX 6.5X10 6 per second. This current distribution represents 
the three major Fourier components of the current profile shown in the inset of 
figure ll(a) . This profile is the same as that in the ion-cyclotron resonance 
apparatus presently being used at the NASA. Lewis Research Center. The principal 
wave number for this distribution is 2koi = 0.140 per centimeter, where 2koi 
is approximately equal to the value of k^ which gives a maximum in the R/k2 
against curve (fig. 5(b)). It is now interesting to examine the contribu- 

tions to relative power absorption caused by the presence of the other Fourier 
components in the current distribution. It can be seen in figure 5(a) 

(R against k^) that when k^ = 4kQ-]_ R is larger than when k^ = 2k 0 j_ (the 
principal wave number). Power transfer, however, is also proportional to the 
magnitude of the current squared; therefore, the contribution to power transfer 
from the kQ^ = 4kop Fourier component must be multiplied by (0.151/0.422 ) 2 
= 0.128. This reduces the effect of the short wavelength component. The con- 
tribution from the long wavelength component (k.Qj = kQ]_) is also small because 
R for , k n = kQi is smaller than R for k n = 2k The current weighting 
factor (0.270/0. 422 ) 2 = 0.630 further reduces the contribution from the long 
wavelength component. Figure ll(a) shows relative power absorption when the 
effects of all three major Fourier components are included. It is seen that 
this curve is not drastically different from the curve in figure 4 where there 
was only one Fourier component. 

Another current distribution which may be examined is the rectangular wave 
distribution shown in the inset of figure 11(b): 

z < a ( 12a ) 

(12b) 

The principal Fourier wave number components for this distribution are 
koi = 0.140 per centimeter and 3kQi = 0.420 per centimeter. Here again, R 
for k n = koi is less than R for k n = 3kQi, but the current weighting factor 
(0.54/l) 2 = 0.292 reduces the contribution from the short wavelength component. 
Contributions from other wavelength components are completely negligible for 
this rectangular wave distribution. Figure 12(b) shows relative power absorp- 
tion for this distribution. 


j* = j( e ik 01 z + 0.54 e 13k 01 z + . . .^ 


J *-0 


z > a 


CONCLUDING REMARKS 

Computer calculations based on the Stix equations (ref. l) for power trans- 
fer to ion-cyclotron waves in a cold collisionless uniform plasma immersed in 
an axial magnetic field and excited by a radiofrequency coil show that radio- 
frequency power transfer bv the plasma waves is relatively independent of den- 
sity (for n e < 5X10 12 cm"*) and increases approximately as the fourth power of 
the plasma radius . It has also been found that for any set of the parameters 
p> s, and n e an optimum radiofrequency coil wavelength exists for maximum 
power transfer. The effect of displacement current for hydrogen plasmas with 
B 0 < 100,000 gauss has been shown to be very small. Also, the lowest order mode 
of oscillation of the plasma was found to contribute about 90 percent of the 
power transferred to the ion-cyclotron waves and to determine the basic shape of 
all the curves. Finally, it should be pointed out that the applicability of the 
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results presented in this report to a particular laboratory plasma depends on 
how well the laboratory plasma is represented by the model upon which the cal- 
culations were based. 


Lewis Research Center, 

National Aeronautics and Space 
Cleveland, Ohio, December 


Administration, 
14 , 1965 . 
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APIEUDIX - COMPARISON OF FOURIER INTEGRAL THEORY WITH 
WITH FOURIER -BESSEL THEORY 

Reference 2 gives the details of a technique employed by Stix prior to that 
of reference 1 for calculating power absorbed by the plasma model discussed in 
this report. To establish the difference between the two methods, calculations 
have been made for the following plasma using both techniques: 

n = 5X10-^* cm" 3 
e 

s = 7.5 cm 



Stix*s equations (ref. 2) contain absolute value signs on P*, but the deriva- 
tion of equation (13a) has been followed through carefully, and there appears 
to be no mathematical reason for inclusion of the absolute value signs. 

If the absolute signs are omitted in equation (13a), there is a value of ^crit 
(see fig. 12) above which the power absorbed by the plasma from the coil is 
negative (ftcrit d s determined by setting the denominator of equation (9) equal 
to zero with k = Fq^). Thus, there is a breakdown in the usefulness of the 
technique of reference 2 at this point. Below ft C rit the results of the two 
calculations are in agreement at only one value of ft. This is the value of ft 
for which k^ = k q}_. In reference 2, the electric field under the finite 
length coil (for 0 < r < p) was represented as a Fourier-Bessel series of non- 
orthogonal (but nearly orthogonal) functions. If' the functions are assumed 
orthogonal in this range, the coefficients of the series can be calculated in 
closed form. This is the technique used by Stix in deriving equation (13a). 

It was thought that this assumption of orthogonality might lead to peculiar- 
ities in the calculations for relative power absorbed. However, Woollett 
(ref. 3) has correctly calculated the coefficients for the Bessel series and 
has found them to agree with the Stix coefficients to within a few percent. 
Thus, the nonorthogonality does not explain the difference between the refer- 
ence 1 and reference 2 techniques. However, the conditions for the two models 
are slightly different. In reference 2 it was necessary to assume that the 
electric field was zero at the left end of the coil when a right running wave 
was assumed for the current sheet. This assumption eliminates a priori the 
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possibility of a left running wave in the plasma. In reference 1 no assumption 
was necessary on the electric field at the ends of the coil. In this case, the 
calculations show that a left running wave exists in the plasma. This wave is 
explained as a reflected wave caused by a mismatch between the natural modes of 
the plasma and the exciting current sheet. When koi = k n> the magnitude of 
the left running wave from reference 1 is zero. This is exactly the point at 
which the calculation on the two models agree. At kgq = k n , the electric 
field is zero at the left end of the coil in reference 1 as well as in refer- 
ence 2 . 

The technique of reference 1 has two other distinct advantages: First, 

it is not necessary to assume s = p as in reference 2, and second, the 
technique of reference 1 can be easily applied to a plasma with a radial varia- 
tion in density by using numerical integration. The reference 2 technique is 
considerably more difficult to apply to a variable density plasma. 
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(a) Plasma radius, 2. 5 centimeters. 


(b) Plasma radius, 5. 0 centimeters. 
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Figure 5. - Power factor R as function of k n . Radiofrequency coil radius, 10 centimeters. 
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Figure 8. - Relative power transfer at maximum P^ ax as function of electron density n e . Plasma radius, 2.5 
centimeters, radiofrequency coil radius, 7. 5 centimeters. 
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Relative power transfer at maximum, P* ax , statvolt/statamp 






Relative power transfer, 




Figure 12. - Absolute value of relative power M as function of Q from 
both references 1 and 2. Plasma radius, 7. 5 centimeters; radiofre- 
quency coil current propagation constant, 0. 140 per centimeter. 
(Method of ref. 2 gives negative P* for Q>^ r ^. ) 
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